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SUMMARY 

The incompressible Navier-Stokes equations are numerically solved for steady flow 
around an ogive-cylinder (fineness ratio 4.5) at angle of attack. The three-dimensional 
vortical flow is investigated with emphasis on the tip and the near-wake region. . 
implicit, finite-difference computation is performed on the CRAY X-MP computer using 
the method of pseudocompressibility. Comparisons of computational results with result 
of a companion towing tank experiment are presented for two symmetric lees, de flow ' cases 
of moderate angles of attack. The topology of the flow is discussed and conclusions are 
drawn concerning the growth and stability of the primary vortices. 

INTRODUCTION 

The motion of a viscous fluid passing by slender body of revolution is an interesting 
flow to study. The complexity of the vortical flow which can arise belies the simplicity 
the body shape. In addition, the fluid-dynamic elements of this flow are common to tho 

found on the components of many flight vehicles. 

The nature of the vortex wake structure which develops on a slender body of re ol 
is primarily dependent on the angle of attack, tip geometry, fineness ratio of the body, a 

to a lesser degree. Mach and Reynolds number. 

On an ogive-cylinder of overall fineness ratio 4.5, (configuration used ,n P 
study) the events that take place as the angle of attack is increase rom 
be associated with a sequence of angle of attack regimes (suggested by Altai and Perkins 
(ref 1)). At zero angle of attack, the flow is axisymmetne. A small increment in ang 
of attack produces a crossflow which eventually causes a steady, symmetric longitudinal 
vL" To for™ on the leeside of the body when the angle of attack exc^a a ce,^ 
minimum The strength of these vorticies, called primary vortices, continues to increase 
with increasing angle of attack. Evidence of the formation of secondary 

sense to the primary vortices may be found near the surface of the ogive. At a cnt^cal^ 
of attack of approximately 35°, an asymmetric disturbance somewhere in th 
a dramatic change in the orientation of the primary vortices to a *’ ^^Tngle 

state For angles of attack greater than 60°, the flow becomes unsteady and as the angle 
of attack approaches 90' a vortex street forms in a manner similar to that foun 

two-dimensional circular cylinder. , 

Side forces approaching the magnitude of the normal force have been attributed 
the highly asymmetric vortex orientations on the leeside of the body. On bodies of grea e 
fineness ratio an array of up to eleven vortices alternating in sign have been observed. The 
precis by which these additional primary vortices form invoives the format, on, growth 
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and shedding of the secondary vortices along the length of the body. The mechanisms 
which cause the flow to become asymmetric are not currently understood. Several hy- 
potheses have been formulated based on the impulsive flow analogy which relates the 
three-dimensional steady flow past a pointed slender body to the two-dimensional un- 
steady flow around a circular cylinder. So far, the analogy has not been fully proven. 

The objectives of this paper are to demonstrate that a valid computation of the sym- 
metric flow around an ogive-cylinder at angle of attack can be achieved in the laminar flow 
regime, and secondly, to provide further insight into the mechanisms which govern the 
growth and stability of the leeside vortical structure. The ultimate goal of this research 
program is to investigate the asymmetric, turbulent flow on bodies of revolution at angle of 
attack and at moderate Reynolds number using a combined computational/experimental 
approach. Successful treatment of the low-Reynolds-number symmetric flows considered 
in this report is intended to be a stepping-stone to these more ambitious attempts. 

METHOD OF SOLUTION 
Governing Equations 

The flow around a slender body of revolution at moderate to high angles of attack and 
low Mach number is a complicated vortical flow with viscous effects playing a primary 
role in the development of the vortical structure. As a consequence, the solution of the 
incompressible Navier-Stokes equations is required. Flow visualization studies reveal that 
the flow about the ogive is steady for angles of attack less than 60°; thus, a time-accurate 
computation is not necessary. For sufficiently low Reynolds number, the flow around 
the ogive is laminar. In a curvilinear coordinate system, the Navier-Stokes equations in 
mass-averaged variables (including the pseudocompressibility terms) can be written in the 
following form: 


dp f du dv dw \ 

Tt ~ 13 \Tx + ai + a~z) = 0 


where 
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The transformation used is of the form 

t = t 

Z = 

T) = T](x,y,Z,t) 

? = c{x,y,z,t) 


with the Jacobian defined as 


J = 


Zz Zy Zz 

r) x rj y T] z 

$z fz . 


and the contravariant velocities 


( 3 ) 


(4) 


U = & + £*u + £ v i> + £ 2 w 
V = Vt + VzU + VyV + riz w 
W = + f z u + + ? 2 u> 

For an orthogonal grid, the viscous flux terms are defined by 


Ev = j{zl + Zl)im— + (/ 2 IZ .. .terms) 

dD 

F v = j(T)l + nl + + (Ry X ---terms) 

Gv = j(C 2 + ? y + + ( fl «- terms ) 


where 


D - jp,u,v.ur J 


( 5 ) 


( 6 ) 
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and 


0 0 0 0 
0 10 0 
0 0 10 
0 0 0 1 

In the above equations, Q is a pseudocompressibility parameter which alters the char- 
acter of the finite-differenced equations to one of a hyperbolic form. This character change 
allows the use of an implicit scheme and relaxes the strict incompressiblity requirement 
during the convergence process. Upon reaching steady state, the effect of/? on the solution 
is negligible. Hence, the steady solution that is obtained is essentially a solution of the 
incompressible Navier-Stokes equations. 

In this investigation the INS3D computer code developed by Kwak, Chang, Shanks 
and Chakravarthy (ref. 2) is applied to solve the equations for flow around an an ogive- 
cylinder body. The numerical approach used in the INS3D computer code is described in 
great detail in reference 2. Briefly, the governing equations including the pseudocompress- 
lbihty terms are written in finite-difference form using central differences for the spatial 
terms. These equations are approximately factored using the factorization technique of 
Beam and Warming (ref. 3). Implicit second-order and explicit fourth-order numerical 
dissipation is added to damp out high-frequency oscillations. The solution is advanced in 
time until a steady state is reached. 

The configuration of the ogive-cylinder model is similar to, and in certain cases, iden- 
tical to configurations used in past experimental studies (ref. 4). The body is a tangent 
ogive-cvlinder of fineness ratio 3.5 with a cylindrical afterbody of fineness ratio 1. The 
5’ 77, ^ curvilinear coordinate system corresponds to the normal, axial, and tangential trans- 
ormed coordinates of the ogive. The tip of the ogive is the location of the origin of the 
(x,y,z) physical coordinates. 

In order to accurately capture the variety of length scales of the flow careful attention 
must be paid to the computational mesh. Near the tip of the ogive, a region of high 
gradients exists owing to the presence of a stagnation point and in addition, because the 
tip is the origin of the leeside vortices. Experiments (ref. 5) have demonstrated that very 
small changes in the tip geometry can have a dramatic effect on the orientation of the 
eeside vortices. In this study, points were clustered at the tip and rear of the afterbody 
using a cosine point distribution. A two-dimensional poisson grid generator (GRAPE, 
ref. 6) was used to distribute the mesh points in a plane passing through the body axis 
ot symmetry. This plane was then rotated about the body axis of symmetry to form a 
three-dimensional wrapped mesh. The grid outer boundary is spherical, and is located 
body lengths from the tip of the ogive. This wrapped grid approach was chosen over 
nree-dimensional grid generation schemes to ensure grid symmetry and also orthogonality. 

A minor drawback of the wrapped grid approach is that the axis around which the grid 
is wrapped becomes a numerical singular line. On this line a one-to-one correspondence 
etween the physical and computational domains no longer exists. To maintain stability 
and accuracy, no grid points were placed on this axis. Figure 1 shows the mesh which was 
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used for this study. The grid is comprised of 49 by 39 by 68 points in the normal, axial, 
and tangential directions, respectively. 

Near the surface of the ogive, the mesh spacing was chosen to resolve the viscous 
effects present in the boundary layers. After extensive study of the effects of grid spacing 
on accuracy, a minimum of seven grid points were placed in the attached boundary layer 
normal to the surface to capture the large gradients present. 

There are several possible sources of error in the numerical method which may affect 
the accuracy of the solution. These errors include the effects of the second- and fourth-order 
artificial viscosity, truncation error, numerical bias, nonresolution of subgrid length scales, 
and boundary condition errors. The effects of the artificial viscosity and truncation errors 
are discussed in references 7 and 8. These effects were shown by analysis and example to 
be minimal, provided that the level of explicit artificial viscosity acting on the pressure is 
allowed to decay during the solution process. To test the resolving capability of the mesh, 
the number of gridpoints was cut in half, and the impact on the solution was observed. 
The effect was found to be minor. The bias errors are thought to be minimal owing to the 
high degree of symmetry in the solution. Estimates of the effects of boundary condition 
errors are difficult to quantify and are a continuing area of research. 

Boundary and Initial Conditions 

Boundary conditions must be prescribed on the six faces of the computational cube, 
and in addition, these conditions must be well-posed and lead to a stable numerical solution. 
The numerical boundary conditions consist of specifying conditions on the pressure and 
velocities on the body surface, at the inflow and outflow boundaries, and near the grid axis 
of revolution. In addition to these four boundary conditions, a circumferential periodicity 
condition is imposed. These conditions, which do not enforce symmetry of the solution, 
are explained in greater detail below. 

Body Surface: 


Inflow: 


u = v = w — 0, 



( 8 ) 


u — U 0 c cos a, v = 0, w = Uoc sin a, p = p x 


( 9 ) 


Outflow: 


Q n i (fman Vr ?) = { <]( f max - 1 j >7? ?) f R 9(£max-2i f] j C )] } (19) 

where q is the vector of dependent variables p, u, v, and w, and R is the extrapolation 
weight factor. 


The extrapolated outflow pressure is additionally modified by the parameter S. This 
parameter is determined by solving a simplified momentum equation and is required to 
insure that mass and momentum are conserved during the extrapolation of the outflow 
velocities. The reader is referred to reference 7 for a deeper explanation of this procedure. 


P {£maxi ■> $) — S p + (£max»hi?) 


where 


S = 


J eii( (momentum corrected pressure)da 
J eiit (extrapolated pressure)da 


( 11 ) 


(12) 


Near the grid axis of revolution: 


Q — {?( £ ? * 7 1 fmai — 1 ) “ 1 “ R [?(£ * ?max) ?(£ ? V ? fmax — 2)] } (^) 


Initial Conditions: 


9°tf,y,c) = 1-0 (14) 

At the beginning of the convergence process, the dependent variables are initialized to 1.0 
in the interior of the computational domain. 

TOWING TANK STUDY 

A companion towing tank study was performed in a 55 gallon tank apparatus (fig. 
2) using a 4-in. long aluminum ogive-cylinder model. The Reynolds number and model 
configuration are identical to those in the computation aside from a small sting which 
was used to support the model. Flow in the crossflow plane was visualized by taking a 
time exposure of the induced motion of a saltwater suspension of neutrally buoyant pliolite 
particles (approximate diam of 0.0035 in.). To improve the laser sheet crossflow plane 
visualization, a cylindrical lens in combination with a short focal-length lens was used to 
produce a 3/10-in.-thick laser light sheet. The model was moved through the water by a 
traversing mechanism and short time-exposure photographs of the induced particle motion 
were taken. 

The particle streaks recorded on the photographs correspond to crossflow instanta- 
neous streamline patterns. The camera and laser sheet were moved with the model to 
produce the instantaneous streamline pattern as would be seen by an observer in the 
model’s frame of reference. After transient effects introduced at the startup of the motion 
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have diminished and the flow has reached steady state, the particle streaks indicate the 
steady streamline pattern. Information concerning the velocity direction and magnitud 
can be obtained from the image by measuring the length and orientation of the particle 

streaks. 


NUMERICAL AND EXPERIMENTAL RESULTS 

Analysis of three-dimensional results is a complicated task owing to the inherent 
difficulties associated with presenting and understanding physical quantnien which . vary 
in three spatial directions. Two different approaches to presenting the data are used l 
this section to cope with the difficulties. In the first approach, data are presented in two- 
dimensional crossflow planes at different locations along the body. Although this approach 
has limitations it is very useful for a flow that is conical or where the flow in the th rd 
dimension is slowly evolving. The second approach involves the three-dimensiona display 
oJTata such as particle path lines. Three dimensional displays of data can sometimes be 
deceptive owing to perspective difficulties, but these displays overcome the limitation of 
the first approach. In the following discussion, some common topological terms (sad e 
node etc. ) are used to describe the flow in the crossflow plane as has been the practice 

for some time (ref. 9). . , , , , 

The two cases computed were for an LjD = 3.5 ogive wit a cy in nca * 
o! LID = 1.0 at angles of attack of 30“ and 45”. The Reynolds number was 1000 based 
on the freestream velocity and overall body length. A comparison of thc comp ““ d 
flow velocity vectors with the corresponding instantaneous streamlines obtained from the 
"owing-tank study is presented in figure 3 (a-e) for a = 30” and l in figure 4 (a-e) for a - 45 
5 denote crossflow-plane coordinates). These figures show the development of the 
flow in crossflow planes at different axial stations along the body. In general, the flow 
development is similar for the two angles of attack. In both cases, the vortices assume a 
symmetric orientation. At X/L = 0.2 (figs. 3a, 4a) no recirculating flo ” rc «'°" ls 
At this station the flow resembles the initial flow observed on an impels ivel Ij s. 
rHmpTVnional circular cylinder. The second station at X/L = 0.4 (figs. 3b, 4b) shows tn 
beginning signs of a recirculation zone, and for the a = 45 c ’ case (fig. 4bj, evi e ^ e ° * 
saddle poin^has appeared on the leeward axis of symmetry. A closeup view of this fl 
is shown in figure 5. In addition to the saddle point on the leeward side, the flow on 
he windward side of the body shows a aero azimuthal velocity po.nt, deveffipment of the 
^flow boundary ,ayer to either side of this point, and ‘ ^£“£ 1 ^ 

it Dasses around the body. The third frame ( X/L = 0.6) of the o; - 30 (fig- ) 
shLs distinctive crossflow separation points, but there is still no evidence of the leeward 
saddle 'formation for this case. At X/L = 0.8 (figs. 3d 4d) the flow "istmcPve 
separation points, saddle points and vortex core centers for both cases. The final fram 
f X/L ~ 0 95) (figs. 3e, 4e), vividly show the primary vortices in their near fully develop 
fa(e Tn addition, the crossflow separation points on the surface have moved in the wind- 
ward' direction creating pockets of nearly stagnant fluid near these points, but no evidence 
of secondary vortices is present in either the experimental or the computational results. 


In general terms, the sequence of events which occurs during the development of the 
vortex flowfields is the same for the two cases studied. The primary difference between the 
cases is the rate at which the vortices develop along the length of the ogive. This result is 
consistent with the impulsive flow analogy which suggests that it should take a longer axial 
distance for the primary vortices to develop on the same body at lower angles of attack. 

The symmetry plane flow (y = 0) over the aft portion of the ogive-cylinder is shown 
in figure 6 for the a = 30° case. Accuracy of this base-flow result may be influenced by a 
small degree of unsteadiness which is not captured by INS3D, but for completeness, this 
result has been included. As expected, a region of 1ow t velocity forms immediately to the 
rear of the cylinder. Fluid originating on the leeside is drawn into this low velocity region 
and is turned towards the streamwise direction by the high velocity fluid which separates 
from the windward side of the cylinder. 

Computed surface pressure plots are shown in figures 7 and 8 for the two cases inves- 
tigated. Considering the a = 45° case (fig. 8). one can deduce several interesting features 
about the vortical flow from the crossflow pressure distribution. As expected, at each 
value of X/L the maximum pressure point is located on the windward ray (0 = 0°). The 
maxima decrease in magnitude progressively from the tip to the end of the body along 
the windward ray. Near the tip of the ogive [X/L = 0.2), the circumferential pressure 
variation is the largest of the five stations considered, and the circumferential pressure 
gradient remains favorable until approximately 0 = 115°. At X/L = 0.95, the circumfer- 
ential pressure variation is the smallest and the onset of adverse pressure gradient occurs 
at approximately 6 - 80 n . The behavior of the pressure at the intermediate stations falls 
between these two extremes. These same trends apply to the a = 30° case (fig. 7) and are 
the expected trends based on the impulsive flow analogy. 

Vorticity contours of the flow field provide further insight into the way in which the flow 
is determined. The axial vorticity contours (w x ) for the two cases computed (at X/L = 0.2 
and X/L = 0.95), are presented in figures 9 (a.b) and 10 (a,b). The intermediate stations 
have the same basic character of these two planes. At the windward ray, the vorticity is a 
minimum and then progressively increases to a maximum and then decreases back to zero 
at the leeward ray. The vorticity maxima are located at approximately 45° on either side 
of the windward ray. These vorticity plots are very reminiscent of typical two-dimensional 
cylinder vorticity contours. 

Three-dimensional representations of the flow field are presented in figures 11 and 
12. Particles were released from points in the windward boundary layer at several axial 
positions. The particles then pass around the body, are entrained into the leeside vortical 
flow, and then they continue downstream. The differences between the two cases are most 
evident in these particle path lines. The vortical flow appears to be stronger for the a = 45° 
(fig. 12) case and also, the vortical flow develops over a shorter distance. 

The corresponding surface skin-friction lines are presented in figures 13 and 14. These 
skin friction lines are the computational analog to surface oil-flow visualizations and were 
produced by releasing particles on the surface and then by following their trajectories. 
Blue particles were released near the windward ray and red particles were released on the 
leeward ray. Magenta particles were used to define ambiguous areas. The red and blue 
particles flow towards each other and form a line where the azimuthal component of the 
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, . f . .. • Thls line which does not emanate from a saddle point, (unlike a 

results presented in recent studies of pointed bodies of revolution (ref. ). 

pnMrT TTQTONS 


■ The flow over a pointed body of revolution ha, been studied from an ' 

understanding the qualitative nature of the flow. 
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Figure 1.- Perspective view of the wrapped computational mesh. 
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MODEL 



Figure 2.- 


Perspective view of the towing tank with model. 
































Figure 4e.- 


Comparison of computed and visualized crossflow velocity vectors 
at X/L = 0.95, q = 45°. 
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Figure 6.- Aft symmetry plane flow at a = 30°. 
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Figure 7.- Circumferential surface pressure in the crossflow plane at 
various stations along the body at a = 30°. 
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Figure 10b.- Axial vorticity contours in the crossflow plane 
at X/L = 0.95, a = 45°. 
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Figure 11.- Computed particle path lines for a = 30 n . 



Figure 12.- Computed particle path iines for a = -45* . 
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Computed surface skin friction lines for a — 30°. 



Figure M,* Computed surface skin iriction lines for o. 45 . 
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